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Abstract. The quantum-mechanical consideration of a passage of fast dime- 
soatoms through matter is given. A set of quantum-kinetic equations for the 
density matrix elements describing their internal state evolution is derived. It 

o 



is shown that probabilistic description of internal dynamics of hydrogen-like 
atoms is impossible even at sufficiently low energies because of the "accidental" 



in 

1 Introduction 

> 

o . 

In the last years interest in studying elementary hadronic atoms (hydrogen-like 
atoms formed by two oppositely charged hadrons) has been growing [1]. 

The basic properties of these atoms are governed by the Coulomb interaction 
between their constituents. The strong hadron-hadron interaction causes only shifts 
of energy levels and increase of widths of s-states of hadronic atoms which are related 



in a simple way to the values of strong scattering lengths of hadrons. 

Thus, the study of hadronic atoms provides valuable information about the low- 



On' 

energy hadronic interactions. At present a number of ongoing or planned experi- 
ments are aimed to measure strong energy level shifts and (or) widths of the ground 
states of the atoms A n - P , A K - p , A n - n +, A n ± K ^ with a high precision [2, 3, 4]. 

The A n - P and Ax- p atoms can be produced at rest in the laboratory frame, 
and a well-developed technique of X-ray spectroscopy can be applied for a precise 
measurement of their radiative spectra. 

The study of the atoms A n - n + and A^kt (dimesoatoms) is a much more 
complicate problem because they practically cannot be produced at rest. The only 
available way to produce a reasonable amount of dimesoatoms is the interaction 
of high-energy projectiles with fixed nuclear targets. The dimesoatoms produced 
in this way move with high (relativistic, as a rule) velocities which determines the 
choice of methods for investigation of their properties. 

For example, the method for measurement of tq (the pionium ground state life- 
time) used in the DIRAC experiment at CERN [4] is based on comparison of the 
observed rate of 7r + 7r - -pairs from the break-up of pionium atoms in the Coulomb 
field of target atoms with the theoretical calculations of this rate which include the 
To as a parameter. 

It is clear that the error in the value of To resulting from such indirect measure- 
ment includes not only experimental errors but also uncertainties of the approx- 



1 Current address: Institut fur Theoretische Physik, Universitat Heidelberg, Max-Plank-Institut 
fur Kernphysik, D - 69029 Heidelberg, Deutschland; e-mail: 01ga.Voskresenskaja@mpi-hd.mpg.de 



1 



imations used in the theoretical description of the pionium production and their 
subsequent passage through the target matter. 

Thus, any improvement of the theoretical treatment of the problems related to 
the DIRAC experiment helps to reduce the resulting error of r measurement. Some 
important results in this direction have been obtained in [5, 6, 7, 8] dealing with 
the theoretical description of the production of dimesoatoms and their interactions 
with separate target atoms. 

In this paper, we will touch on a problem of describing the internal dynamics of 
dimesoatoms moving through matter and interacting with its atoms. It is commonly 
supposed [9] that the dimesoatomic dynamics can be described with the help of a set 
of kinetic equations for the probabilities Pi(z) to find the dimesotom in the definite 
quantum state \i) at the distance z from the production point. 

Although the form of these equations is similar to the Pauli quantum kinetic 
equations [10], we will refer to them as the classical (probabilistic) ones because 
they ignore interference between different states of dimesoatoms during their passage 
through the matter which is pure quantum effects. 

It is well known that these interference effects play a decisive role in giving 
rise to the so-called superpenetration phenomena [11] at ultrarelativistic energies. 
However, at the typical dimesoatom energies in DIRAC most interference effects are 
strongly suppressed by the time-formation effects. 

But in the case of dimesoatoms some of interference effects occur even at suffi- 
ciently low energies because of accidental degeneracy of energy levels of hydrogen-like 
atoms [12]. Since these effects cannot be treated in the classical approach of the pa- 
per [9], one have to use a more sophisticated approach for consideration of internal 
dimesoatom dynamics based on the density-matrix description of quantum systems. 

The main goal of this paper is to derive a set of equations for the density ma- 
trix elements describing the internal dynamics of multilevel atoms moving through 
matter and interacting with its atoms. Although for simplicity the derivation is 
restricted to the case of elementary (two-particle) atoms, the method can be easily 
generalized to a more complicated case of arbitrary composed system. 

The plan of the paper is as follows. In Section 2, the main approximations are 
described. In Section 3, a set of equations for density matrix elements is derived. 
In Section 4, the correspondence between density-matrix approach and probabilistic 
approach of the paper [9] is discussed. In Section 5, a brief summary of the main 
results is given. 

2 Main approximations 

Any discussion of relativistic dimesotom interactions faces, first of all, the problem 
of the relativistic description of a bound system. Application of usual methods based 
on the Bethe-Salpeter equation [13] would lead to overcomplication of the technical 
part of the problem under consideration, which may mask the main goal of this 
work. 

To avoid this difficulty, we will use the following trick. We will first consider 
the interaction of fast but still nonrelativistic dimesoatoms with the target. In 
this case the relatively simple mathematical formalism of nonrelativistic quantum 
mechanics can be used for the description of both internal and external dynamics 
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of dimesoatoms. This results in simple energy dependence of the density matrix. 

The next step is the "analytical" extension of this nonrelativistic result to rela- 
tivists energies appealing to the physical intuition. Although this way of derivation 
cannot be considered as rigorous we believe that it allows one to get the most es- 
sential features of the desired result. 

We restrict our consideration to the case of rather high dimesoatom velocities 
and rather thin targets so that the "frozen" target approximation can be used. In 
other words, we will suppose that the positions of target atoms practically do not 
change during the passage of a dimesoatom trough the target. 

In this case, the stationary Schrodinger equation can be used to describe the 
interaction between the dimesoatom and the target: 

H^ = Etjj, (1) 

H = H int + H ext , (2) 

H int = -^-A r +V int (r), (3) 



H ext = -^A^ + V ext (f, R; {p}) . (4) 



1 
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Here = ?/>(r + ,r_), r + and r_ are the radius- vectors of the positively and nega- 
tively charged mesons respectively, m + and m„ are their masses; M = m + + 
\i = m + m_/M; R = {m + r + + m_r_)/M is the radius-vector of the dimesoatom 
center-of-mass; r = (r + — r_) is the radius- vector of relative motion of dime- 
soatom constituents; {p} is the set of the radius- vectors of the target atom center-of- 
mass; Vi n t(r) is the potential energy of interaction between dimesoatom constituents; 
V ext (f, R; {p}) is the potential energy of interaction between the dimesoatom and 
the target atom: 

V ext = e5>(f + - p k ) - W- ~ Pk)] , (5) 
k 

r + = R + ^r, f. _ = R — rjf , (6) 
, rn- m + 

Z = W , V = W = 1-Z, (7) 

$ is the electrostatic potential of the target atom at rest; e is the elementary charge. 

Let us consider that the beam of dimesoatoms is in a definite state of the discrete 
spectrum \i) and moves in the positive direction of the Oz-axis (p = (0,0, p) is the 
beam momentum). The target is supposed to be a plate of thickness L and of the 
infinite transverse size, placed perpendicular to the beam, so that one surface has 
z = and the other has z = L. Then, it the z < 0, 

i>(r+, r_) = exp(pR) (8) 

H int (j)i(r) = Ei<pi(f) , (9) 

where £j is the energy of the state 

The solution of (1) in the region z > will be sought in the form 

^(f+, f_) = exp(ipR) ■ F(R, r, {p}) . (10) 
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Then the equation for F following from (1) is conveniently written as 

= (-£jfc + H M + vJ)F(z,b,r>,{Pl), (11) 

— * — * 

where z and b are the longitudinal and transverse parts of the radius-vector R = 
{b,z}. 

At rather high dimesoatomic energies, such that the condition 

pa>l, a ~ (m e • a ■ Z 1 ^)' 1 (12) 

is satisfied, the term d 2 F/dz 2 on the left side of (11) can be neglected [12] (in (12) 
a is the screening radius of the target atom, Z is its atomic number and m e is the 
electron mass). 

In this approximation the quantity 

F = exp (^- £ izj ■ F ( 13 ) 

obeys the equation 

dF 

-i^ = H-F, (14) 

1 M 1 

H = -— A s + —H mt + -V ext (15) 
2p p v 

which looks like the nonstationary Schrodinger equation. 

In (15) v is the dimesoatom velocity, which in the nonrelativistic limit coinsides 
with the ratio p/M. But keeping in mind the subsequent "analytical extension" 
of the below results to the region of relativistic energies, we will assume not the 
equality p/M = v in the general case but the relativistic relationship p = Mvy. 
With this prescription the right energy dependence of the time-formation effects is 
reproduced (see also [14]). 

3 Density matrix kinetic equation 

Let us write the expression for the density matrix describing the dimesoatom state 
after covering the distance z in the matter 

p(z; n, f 2 ) = S- 1 J dbdb[db' 2 dr[dr' 2 (G(z; b, f ± ; b[,f{; {p}) 



xG*(z; b, r 2 ; 6£, {p}))^ 0^)0*^') . (16) 

Here G is the Green function of (14), the solution of which is given by the expression 

F(z, b, f, {p}) = j G(z; b, f; b', f '; {pj^r'W'db' , (17) 

E is the beam cross section area, and the symbol ()^ means averaging over all 
possible states of target atoms. 
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To perform this averaging in the explicit form, we will use the path integral 
representation [15] of the Green function: 



G(z; b, r; b, f ; {p }) = J Db(z)Dr(z) exp(iS) , 
where the action S of the dimesoatom interacting with the matter is 



(18) 



p ( db(z') 
2 I dz' 



wyp ( dr(z') 
~Y \ dz' 



o 

-— V int (f(z')) - -V ext (V; b(z'), f(z')) 



(19) 



The path integration in (18) is performed over all trajectories |&(z), f(z)| beginning 
with the point {&(0),f(0)} = jft', r'j and ending with the point {b(z),r(z)} = 
{b,r}. 

It was shown in [16] that 
r 1 

J dz'-V ext (z';b(z'),f(z')) = *E&(z-z k ){ X (b(z k )+S3{z k )-T k ) 

k 

- x (b(z k )- V s(z k )-r k )} , (20) 

where the Heavyside step-function -&{z) is for z < and 1 for z > 0; s and f k are 
the transverse components of the vectors r and p k \ z k is the longitudinal component 
of the vector p k . 
The quantity 



— oo * ' 

Bz = b + £s-r k , B v = b-rjs-T k , 



(21) 



is the phase shift acquired by the wave function of the hadron while it moves through 
the field of the isolated target atom. 
Let us use the relation [11] 



cxp 



L k 



{p} 



exp |-n J dp(l - f(p, {*}))} , (22) 



where no is the number of target atoms in the unite volume; f(p k , {x}) is the arbi- 
trary function of the target atom coordinates and of the set {x} of others variables. 

Then the average of the product of two Green functions over the target atom 
coordinates is 



G(z-b u ri,b[,r[-{fi)G*(z-b 2 ,r 2 -b^r^{p}) 



{p} 



= J Db^Dhi^Dn^Dniz) exp (iSi - iS 2 - n J dz' J cudA , (23) 
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Sip) = / dz' 





, (24) 
(25) 



Axi( 2 ) = X pi(2) (* ') + f si( 2 ) (2 ') - r] - x pi(2) (2 ') - VSi(2) (z')-t]. (26) 
With the replacement 



r^r + b + (z'), 

where the f is impact vector and 

b + (z') = (b 1 (z') + b 2 (z')-f) , 
it is easy to see that u depends only on 

b_(z') = (b 1 (z')-b 2 (z')-r) 



(27) 
(28) 

(29) 



and does not depend on b+(z'). 

This allows to perform path integration over the variables bi(z) and b 2 (z) i n (23) 
in the explicit form (see e.g. [16]) and to represent the density matrix (16) in the 
form 



p(z;fi,f 2 ) = J Dfi(z)Df 2 (z)df[df 2 



where 



x exp { iSi - iS 2 — n J dz'Vt) 0i(r/)0*(r 2 ') , 



Q = J codr = n(s 1 (z'),s 2 (z')) , 
s 2 ) = Jdr {r(f, sO + r (r, s 2 ) - T(r, sjT* (f, s 2 )} , 

r(f,si( 2 )) = 1 -exp [ix{T - £si( 2 )) ~ ix{? + VSi(2))} ■ 



(30) 



(31) 



(32) 



(33) 



r(r, s) is the interaction operator of the Glauber theory for interaction of dime- 
soatoms with target atoms [17]. In particular, the transition amplitudes of the 
dimesoatoms between the states \i) and \k) in the Coulomb field of target atoms are 
related to T(r, s) by the equation 



A ki (q) 



2tt 



rf 2 rrf 3 re^*(r)^(r)r(f,s) 



(34) 



Path integrations (like the usual ones) in (30) can be performed by numerical 
methods ([18]). 

But for the purpose of the qualitative analysis it is more useful to deal with the 
kinetic equation for the density matrix 



.dp{z;fi,f 2 ) 
dz 



i>7 



[Hintin) - H* nt (r 2 )\ p(z; r u f 2 ) 



-m fi(si, s 2 )p(z;fi,f 2 ) , 



(35) 
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which is the consequence of the path integral representation (30) (see, for example, 
[15]). 

Because of the hermiticity of Hi nt this equation can be rewritten in the operator 
form 

*7T" = ~~ \. H i"*> P\ ~~ in o^P > ( 36 ) 

oz vy 

where the first term on the right side (Liouville term) describes the causal part of 
the internal dimesoatom dynamics. The nature of the second term is pure stochastic 
because of the stochastic distribution of the atoms in the target. 2 
In the dimesoatom rest frame Eq. (35) takes the form [14] 

.0p(*;fi,f 2 ) 



dt 



= H int (fi)p(t; fi, f 2 ) - H* nt (f 2 )p(t; fi, f 2 ) 

-ivynoQ(si,s 2 )p(t;fi,f 2 ) (37) 



which is similar to the form of Eq. (116) of paper [19] for the density matrix of 
atoms moving in the laser fields. But the meaning of the terms on the right side of 
this equation and on the right side of our Eq. (37) is different. 

The causal part of (116) also includes the term describing interaction of atoms 
with the laser fields while the stochastic term of this equation describes effects of 
spontaneous relaxation caused by interaction of atoms with the quantized electro- 
magnetic field. 

Strictly speaking, a similar term should be added to the right part of (37) too. 
But it can be shown that its influence on the internal dynamics of the dimesoatoms 
is negligible compared with similar effects caused by interaction of dimesoatoms 
with target atoms. 3 For this reason we do not include the effects of spontaneous 
electromagnetic relaxation in our consideration. 

There is effect that play important role in the dimesoatom dynamics which is 
not taken into account in (37). It is the effect of instability of dimesoatoms caused 
by the possibility of their annihilation into neutral hadrons. We will include it 
phenomenologically. 

With this aim let us consider the density-matrix elements 

p lk {z) = J $W*(f>(*; r, r')drdr' . (38) 
Taking into account the Schrodinger equation 

H int 1pi = Si^i (39) 

it is easy to see that matrix Eq. (35) is equivalent to the following system of equations 
for the density-matrix elements: 

= i\kPik ^ik,lmPlm , (40) 

l,m 



2 For each separate projectile the target "is frozen" in the definite state. But because of the 
stochastic distribution of the projectiles in the beam these definite "frozen" states of the target 
are distributed stochastically too. 

3 Interaction of dimesoatoms with the quantized electromagnetic field causes only deexcitation 
of excited dimesoatoms through emission of real photons. Interaction of dimesoatoms with target 
atoms causes both excitation and deexcitation of dimesoatoms through exchange of virtual photons 
with target atoms. 
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where 



A ife = (e k - £i)/u7, 



(41) 



It is evident that after including the instability effects the solution of this system 
at n = ("empty" target) 



Pik{Z) 



p ik (0) exp[iA ik ■ z\ 



is modified to 



Pik(z) = Pifc(0)exp 



t>7 



i(e k - Si) - -(Ti + T k ) 



(42) 



(43) 



where are the widths of the states \i(k)). 
Hence Eq. (40) must be replaced by 



dpi k (z) 
dz 



1 

i>7 



i(£ k -£i) - -(Ti + T k ) 



pik (z) -n J2 ^ikMPim(z) . (44) 

l,m 



For the subsequent discussion it is convenient to represent (44) in the integral 
form 



Pik = Pik(0) exp[-a ife ■ z]-n J2 ^ik,im J exp[-a ik ■ (z - z')]pi m (z')dz' , (45) 

(46) 



l,m 



&ik — 

t>7 



i(e k - si) - -{Ti + r fc ) 



4 Time-formation effects 

It is easy to check that the coefficients Vt ik ^ m (41) are different from zero if the 
relations 

mi~m k -mi + m m = , h - l k - h + l m = 2s (47) 

are satisfied. Above, m,i(k,i,m) an d h(k,i,m) are the magnetic and orbital quantum num- 
bers of the states \i(k,l,m)) and s is an arbitrary integral number (the quntization 
axis is chosen to coinside with the Oz-axis). 

Since the dimesoatoms are produced practically only in the ns-states (I — 0, m — 
0) [4], the "selection rules" for the quantities Pi k (z) 



Pik(z) 7^ if and only if rrii = m k , Zj = l k + 2s 



(48) 



follow from (42), (47). 

No restrictions on the values n^ k ) — principal quantum numbers of the states 
\i(k)) — are imposed by the dynamics of interaction of dimesoatoms with target 
atoms represented by the quantities Qik,i m - But they are imposed by the so-called 
"length (time)-formation effects", 4 represented by the oscillating factors exp[iA ik ■ 
{z-z')\ in (45). 

4 These effects belong to the wide class of the effects systematic investigation of which began in 
the papers [20, 21] (further development see in [22]). 
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At super-high energies, such that A^L <C 1 (L is the thickness of the target), 
the factors exp[iAik ■ (z — z')\ ~ 1 for any values of the indices i, k, and the relations 
between the values of the diagonal and off-diagonal density matrix elements Pik{z) 
are governed only by the properties of the coefficients flik,i m - 

For "low" energies (such that A ik \ 3> 1, where A is the typical dimesoatom 
free-path length in matter) the numerical values of the off-diagonal density matrix 
elements Pik{z) (i ^ k) become negligible compared to the diagonal ones if Si ^ Sk- 

If energy levels of the system under consideration are non-degenerate (ei ^ £k, if 
i ^ k), only the diagonal density matrix elements pa = Pi survive at "low" energies 
because of length-formation effects, and the internal dynamics of such systems can 
be described by a system of kinetic equations for probabilities Pi(z) 



which coinside with the system equations of [9]. 

Taking into account the "selection rules" (48), we clearly see that this statement 
is also valid when there is only degeneracy with respect to the value of the mag- 
netic quantum number that is, if the hadron-hadron interaction has only the central 
symmetry. 

However, the Coulomb interaction in the hydrogen-like atoms has additional 
symmetry that causes extra (accidental) degeneracy of the energy levels with respect 
to the value of the orbital quantum number. Because of this, even at very low 
energies some off-diagonal density matrix elements Pik(z) with 



survive together with diagonal ones. This excludes the possibility of describing 
of internal dynamics of these atoms in the framework of "classical" (probabilistic) 
approach [9] at any energies. 5 

5 Conclusion 

The accurate description of internal dynamics of relativistic dimesoatoms moving 
through matter requires the density matrix formalism. The kinetic equation for the 
density matrix could be represented in the form of the Liouville equation only if the 
positions of the target atoms will change in time according to a definite law. Because 
of the stochastic change of these positions, the part of the kinetic equation for 
density matrix describing the interaction of dimesoatoms with target atoms becomes 
more complicate. The main result of the present paper is the derivation in the 
explicit form (32) of the term describing this interaction. That fact that we have 
automatically obtained an expression for (32) in terms of interaction operators of 
the Glauber theory can be considered as additional confirmation of the validity of 
the approximations used. 

5 Strictly speaking, strong interaction between dimesoatomic constituents eliminates the acci- 
dental degeneracy of the dimesoatomic energy levels. But this effects is so small numerically that 
it has no significance for the problem under consideration (for example, (e2 P — £2s)/ £ 2s ~ 10 -3 ). 



dz 



cuPl(z), en = -nofluji - a u 5u 

i,l 



(49) 



— £-ki li — ^fc ~r" 2s 



(50) 
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